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Abstract. The objections raised by Vilar and Rubi [cond-matlarXiv:0707.3802Vl] 
against the definition of the thermodynamical work appearing in Jarzynski's 
equaUty are shown to be misleading and inconsistent. 

Quid dices de primariis huius Gimnasii philosophis, 
qui aspidis pertinacia repleti, licet me ultro dedita 
opera millies offerente, nec Planetas, nec J, nec 
perspicillum videre voluerunt? Verum ut ille aures, 
sic isti oculos, contra veritatis lucem obturarunt. 

Galileo (in a letter to Kepler, 1610) 



1. Introduction 

In a recent post, Vilar and Rubi (VR) Jj ascribe to Imparato and Peliti [2j the claim 
that the standard definition of work 

Work = Force x Displacement, (1) 

should be unconditionally replaced by a 'nonstandard' definition 

dWiv^-xAf, (2) 

in which the force and the displacement have their role interchanged, when considering 
the work performed by a force on a system. They argue that "this 'nonstandard' 
definition of work is thermodynamically inconsistent at both the microscopic and 
macroscopic scales and leads to non-physical results, including free energy changes 
that depend on arbitrary parameters" . The dispute arose from the claim set forth 
by Vilar and Rubi in a previous post [3], in which it was argued that the connection 
between the microscopic work W performed by a time-dependent force on a system 
cannot be used to estimate free energy changes. 

In the present note I shall argue the following points, which are already clear to 
any honest reader of ref . [5] : 

(i) The expression ([2]), surprising as it is, is a straightforward consequence of the 
standard definition of the thermodynamical work performed on a system, for the 
special case considered in [5], namely, when a uniform but time- varying force is 
applied to a particle subject to a given potential; 
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(ii) This expression yields, via a straightforward application of the First Principle 
of thermodynamics, a correct evaluation of the free-energy change of a 
thermodynamical system undergoing a reversible transformation; 

(iii) The 'inconsistencies' claimed by VR to be produced by this expression of the work 
correspond to bona fide energy differences which have observable consequences. 

I shall also argue that VR's confusions stem from the fact that the thermodynamical 
work on a system represents the work done by the system one considers on the 
external bodies which act on it, rather than the work done on the system itself by the 
external bodies: a point stressed, e.g., at the beginning of Gibbs's foimding book [3] 
on Statistical Mechanics, and that VR fail to appreciate. 

I shall first discuss these points in the context of equilibrium thermodynamics. 
Further points are relevant when considering manipulated systems, in particular small 
systems for which fluctuations are important. 



2. Reversible work on the harmonic oscillator 



Let us consider a simple thermodynamical system, i.e., a one-dimensional oscillator 
characterized by its mass m and spring constant k, kept at a fixed temperature T. 
The system is described by the hamiltonian 



Hip,.) = ^ + lk.^. (3) 

In the following we shall focus only on the displacement degree of freedom, namely x. 
Its equilibrium distribution is given by 

p-kx^ /2kBT 
P'-'i^) = z ' 

where Z is given by 

Jdx e-'=-'/2'=BT ^ ^2nk^T/k. (5) 

We shall now apply a uniform, but time- varying, force f{t) to the system. We 
wish to evaluate the thermodynamical work performed on it, as the applied force 
changes from /o = to /, so slowly, that the system can he considered to remain 
at thermodynamical equilibrium at all times. This is called the reversible work in 
thermodynamics. 

Following the method described by J. W. Gibbs [4J and R. C. Tolman [5], one 
proceeds as follows: 

1. One writes down the hamiltonian of the system in the presence of the applied 
force: 

H{x,f ) = hx^ ~ f{x~^). (6) 

Here 7 is defined as the point in which the potential of the applied force vanishes. 
This point might depend on /, but we shall momentarily assume that it is fixed. 
It is determined by the actual device used to apply the constant force on the 
system, as discussed in the following. 
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2. One applies either Gibbs's equation (117) 4, p. 45], or Tolman's equa- 
tion (124.1) [51 p. 542], to obtain the therniodynamical work dW associated with 
a small variation d/ of the applied force: 

dW = (^^^ df = - {{x - 7)) d/ = - ((x) - 7) df. (7) 



{A}^- I dx A{x) e-f^(^./)/teT^ 



In this equation, {A) is the canonical average of the function A{x) 
1 
Z 

In our case, one obtains 

(^) = {, (9) 
from which dW can be calculated via equation ([7]). 

One integrates the result with a variable force /' from the initial value /o = to 
the final value /, obtaining 




d/'^-/'({-7)d/' = -g+7/. (10) 



/'/ 

In this expression, AF is the change in the Helmholtz free energy, F ^ E — TS. 
Since it is easy to see that in the present system the entropy S does not change 
during the manipulation, we can equate it with the change in the internal energy 
E. We have therefore 

AE^-Q+^f. (11) 

4. Since the average value of the applied potential is given by 

{U) = -fi{x)-j), 

by subtracting it from the above result, we obtain the change of the energy of 
self-interaction of the spring 

AE'^\f) ^ AE - (U) ^ ^. (12) 

It would not be necessary to consider this elementary exercise in statistical 
mechanics, were it not for the fact that in their recent post J. Vilar and M. Rubi [T] 
(objecting to a similar derivation contained in 2J) have found that this result 
is "inconsistent and unphysical both at the macroscopic and microscopic level." 
Equation ([7]) is the one that VR incriminate. The two authors are chagrined by 
the following facts: 

1. Let us first consider 7 = 0. Then the free-energy change ((TO)) is negative. Now, 
non-spontaneous processes should lead to positive free-energy changes. This is in 
contrast with previous results, including ones on macromolecules [5]. Moreover 
this result holds for any system described by the hamiltonian including a 
macroscopic spring. This is in contrast with the results of elementary physics. 

2. Moreover, VR claim that the parameter 7 does not have any physical 
interpretation, and that therefore in this result the free-energy change does not 
depend on the actual physical system but rather on its mathematical description. 
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VR notice that the free-energy change given by equation pop is negative since it 
also contains the potential energy associated with the external force. They claim 
that it is inconvenient to have the particular properties of the applied external 
force embedded into the results. Therefore they deviate from the definition of the 
thermodynamical work given by Gibbs, who explicitly states [H p. 4, footnote] that 
the energy function of the statistical system should include "that energy which might 
be described as mutual to that system and external bodies" . It lies on them, therefore, 
the burden to show that their 'nonstandard' definition of the thermodynamical work 
is preferable to Gibbs's and Tolman's one. They shun this burden by failing to notice 
it. 

They should however agree that, if the potential energy of the interaction of 
the system with the bodies that provide the constant force is taken into account, 
both objections raised above disappear. VR proceed instead as if the expression (fTTj) 
contained only the energy of interaction of the system with itself. 

I now show how the result pT|) corresponds to the variation in the total energy 
(as defined in the above text by Gibbs) when the external force is applied by two 
physically reasonable devices. I shall then discuss why the apparent paradox of point 
1. is such only in the minds of the authors of ref. [I] and their followers. But I now 
wish to stress the point which probably lies at the heart of VR's confusion, by quoting 
at length from Gibbs's treatise. 

Returning to the case of the canonical distribution, we shall find other 
analogies with thermodynamics systems, if we suppose, as in the preceding 
chapters!! that the potential energy (eg) depends not only upon the 
coordinates qi . . .q„ which determine the configuration of the system, but 
also upon certain coordinates ai, 02, etc. of bodies which we call external, 
meaning by this simply that they are not to be regarded as forming any part 
of the system, although their positions affect the forces which act on the 
system. The forces exerted by the system on these bodie^ will be represented 
by —deq/dai, ~deq/da2, etc., while —deq/dqi...~deq/dqn represent all the 
forces acting upon the bodies of the system, including those which depend 
upon the position of the external bodies, as well as those which depend only 
upon the configuration of the system itself. It will be understood that ep 
depends only upon (71, . . . (7„,pi, . . .p„, in other words, that the kinetic energy 
of the bodies which we call external forms no part of the kinetic energy of 
the system. It follows that we may write 

de dtq ^ (104) 
da\ dai ' 

although a similar equation would not hold for differentiation relative to the 
internal coordinates. 
Thus Gibbs's expression of the elementary reversible work 



(where, in Gibbs's notation, the bar denotes the average over a canonical distribution) 
represents the average work done on the external bodies by the system (with changed 
sign), and therefore, in particular, does not vanish even if the coordinates of the system 
do not change over the time interval considered. 




(13) 



X See especially Chapter I, p. 4 (Note by JWG). 
§ My italics (LP). 
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Electrostatic device 



To illustrate this point, let us set up a device for applying a uniform but time- 
dependent force on our oscillator. We can use, for instance, the following electrostatic 
device. Let us assume that the mass of the oscillator carries a small charge q. We 
take two point-like bodies at infinity, one with the charge +Q and the other with the 
charge —Q. We then let these two charged bodies come closer and closer to the origin 
(the equilibrium point of the oscillator) , by letting the charge +Q be situated at the 
point —X + 7, and the charge —Q at the point X + j. Thus the electric field acting 
on the oscillator at point x is given by 



E ^ 



Q_ 

47reo 

Q 



1 



1 



(x. 



X-7)2 
(a; -7)2 



{x + X 
f X2 



2^60 [{x - -f)^ + X^]^ ~ 
Q f 1 3(x-7)2 
27reo \X^ X^ 



A{x - 7)2^2 

5(x-7)4 
I — 1^ 'J. — I 

X6 



(14) 



If X is large enough, then all terms beyond the first one are negligible, for the expected 
excursions of the oscillator from the origin. Then the force applied by the charge Q is 
given by 

qQ 1 



f = 



(15) 



27reo X2 ■ 

Let us choose Q such that, even for the largest force /i which we wish to apply, X 
is so large that the terms beyond the first in equation (|14p are negligible. Thus by 
moving the charges ±Q from infinity to ±X -I- 7, always symmetrically around the 
point 7, we can apply a uniform but time-varying force to our oscillator. It is now 
clear that 7, far from being a fictitious parameter, corresponds to the location of the 
center of the device by which a uniform force is applied to the system we are studying. 
In order to change 7, external work must be supplied to the apparatus. 

Let us now evaluate the internal energy of the system as a function of X. We 

have 



E : 



kx^ + U{x,X) 



- I dx e-^(-'^)/'=BT 



-kx' 
2 



1 



1 



qQ 

47reo \x + X — J X — x + ^ 



The first term yields 













2 



{ix-{x)r) + {xy 



= -knT 
2 



2 k ' 



(16) 



(17) 



The first term is given by the equipartition theorem, and the second by equation ^ 
One can expand the second term in powers of 1/X, obtaining 

qQ 

' ' X^ 



{U{x,X)) = - 



2n€o 



(18) 



Thus, if {{x - 7)2) /X^ < 1, we have 

qQ ii^-i)) 



{U{x,X)) 



X2 



11 
k 



if, 



(19) 



Comment on: Inconsistency of the nonstandard definition of work 



6 



where we have exploited and (O. Summing up, we obtain 

^^ = ^fcBr-^ + 7/, (20) 

in agreement with equation (jlip . 
Gravity-field device 

A simpler conceptual experiment can be set up imagining that the oscillator mass is 
constrained to move along a line, which can be rotated in the vertical plane. Let m 
be the oscillator mass, g the acceleration of gravity, and let the hinge be placed at 
a; = 7. If the line is now rotated clockwise by an angle 6, the oscillator mass will acted 
upon by a uniform force, directed towards increasing values of x, and of intensity 
mi? sin ^. On the other hand, if the mass is at location x, its height with respect to 
the horizontal line passing through the hinge is given by z = — (a; — 7) sin 9. It is then 
a simple matter to evaluate the average oi U{x,9): 

{U{x,e)) = mg{z) = -mg sine {x--f) = -^+-ff. (21) 

Adding to it the average elastic energy ifc^x^) we recover equation (fTTj) again. But 
it is amusing to verify that this result corresponds indeed to the work done by the 
system on the external device. Let us consider the line to be tilted by 9, and the 
position of the oscillator to be x. Then the oscillator applies to the rectilinear guide 
a torque 

T = mg cos9 {x — j). (22) 

As the angle changes by d9, this torque executes on the guide a work 

T d9 = mg{x — 7) d sin 9. 

The reversible elementary work made by the system on its environment is given by 
the average of this expression, namely 

- AW'"''' = (r) d9 = mg {{x) - 7) d sin 61, (23) 

where, according to equation ([9]), (x) = f /k = mgsai9 /k. The change in the internal 
energy due to the transformation is given by dW^'^^ , integrated between and the 
final value of 9. It is easy to check that it yields again the result (jlip . 

When the rectilinear guide is tilted, the oscillator spring is stretched and its elastic 
energy is increased. On the other hand, the potential energy of the mass in the gravity 
field can either increase or decrease, and the resulting total energy change can be of 
either sign. If 7 = 0, one has, for instance 

VR claim that this result is inconsistent, because a negative free-energy change 
(which coincides in our case with the energy change) would imply that the process is 
spontaneous, and that the spring is unstable, in contradiction with elementary physics. 
They fail to notice, however, that, if the rectilinear guide is free to rotate around the 
origin, the system is indeed unstable: the guide would rotate till it reaches a vertical 
stand, with the oscillator mass hanging on the spring. Thus, far from being unphysical, 
the result yields the correct prediction for the physical setup one is considering. Of 
course, in an actual experiment, one would constrain the guide at a given angle 9, and 
the oscillator will find equilibrium around a point {x) given by equation ([9]). 
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3. Reversible and fluctuating work 

The textbook definition of reversible work is the work performed when the 
thermodynamic transformation is so slow that the system can be considered to stay at 
thermodynamic equilibrium at all times. In this case, the trajectory average coincides 
with the ensemble average, at least if equilibrium statistical mechanics holds. Then 
the performed work does not fluctuate, and one trivially has 

p^rcv ^ ^^rcv^ ^ -fceTlog (^e'^'°^ /''^^'^ . (25) 

VR claim that this equality does not hold, presumably because in their mind the 
reversible work (which is a canonical average) fluctuates. On the other hand, a clear 
distinction was made in ref. [5] between reversible and fluctuating work, a distinction 
that VR chose to ignore. For the benefit of the reader, I recall the definition of 
the infinitesimal fluctuating work on a system whose microscopic state is denoted 
by a; = {xi), and described by the hamiltonian H{x,^), depending on an external 
parameter fj.: 

dW=^^^d,. (26) 

We then have, for a given infinitesimal change d/^, 

dW^--"^ = (dW) , (27) 

where the average is taken with respect to the canonical distribution with the given 
value of /z. Notice that the fluctuating work does not depend on the change in the 
microscopic state x of the system, but on the change of the external parameter /x, 
because it represents the work done by the system on the external bodies that act 
on it. One can thus understand why it does not vanish if /i is suddenly changed: 
if, e.g., we suddenly push the charges ztQ closer to the origin in the electrostatic 
device, we have to provide some work, part of which changes the interaction energy 
of the oscillator with the charges. By the same token, if we change 7, e.g., by rigidly 
displacing the field-creating charges Q, we have to provide work on the system, even 
if the oscillator's mass does not move. 

The distribution of the fluctuating work exhibits a number of interesting 
properties, among which the remarkable equality ^e"'^/'^^'^) = q-^^/'^'bT j^gj-ived 
by Jarzynski [71, and which Hummer and Szabo [H] showed how to exploit in order 
to obtain information on the equilibrium free-energy landscape from nonequilibrium 
experiments. VR object to this development, claiming that the above definition of the 
fiuctuating work is unphysical and inconsistent. We have just seen how nicely it fits 
with equilibrium statistical mechanics, as defined by Gibbs and explained by Tolman. 
However, other quantities also exhibit remarkable distributions. Let us consider a 
system described by the hamiltonian 

H[x,^i) = H^[x)-^lF[x). (28) 
Then the fiuctuating work defined above is given by 

AW = -F{x) d^, (29) 
and satisfies Jarzynski's equality. On the other hand, we can also define the work dWo 

by 

dW^ = ^lY, — dx„ (30) 
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which represents the work done by the environment on the system. As recently 
discussed by Jarzynski |9J in more detail, this work satisfies an identity found long 
ago by Bochkov and Kuzovlev 'lOj, namely 



However, it is true that it is difficult to exploit this identity in order to recover 
information on an equilibrium quantity like AF. Indeed, what VR have brilliantly 
shown in [3] is that Wo cannot be used to reconstruct free-energy landscape, but they 
fail to inform the reader of [T] that their arguments concern the use of Wq, leaving 
the impression that their objections concern W and the use of the Jarzynski equality. 
Now there is no problem in applying the Jarzynski equality to W. The resulting AF 
contains a contribution from the interaction between the system and the environment 
which, contrary to VR's statements, is easily subtracted off (see, e.g., the "histogram 
method" discussed in [HI [12]). It is the responsibility of the researcher to choose the 
most appropriate tools for one's task. One should choose a spoon to eat one's soup 
and a spade to dig a hole: VR appear to prescribe everybody to pick up the spoon 
and then they lament that it is not possible to dig holes. 

4. Conclusions 

We have seen that VR's objections against ref. [5] stem from a biased and misleading 
reading of it, and from their failure to appreciate some basic concepts in statistical 
mechanics. I am at a loss to understand why as serious and competent physicists as VR 
could fall in such blunders, unless their confusions arise from an aprioristic hostility to 
the recent exciting developments in the statistical mechanics of manipulated systems. 
In this case, they would remind of Galileo's colleagues, cited in the letter I have posted 
in limine, who refused to look in the telescope because it did not fit within their world 
view. If it is so, let them be happy to encourage their followers to raise objections 
based on even faultier arguments than their own 13J. I shall have no more to say on 
their subject. 
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